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^ ■ Abstract 

. We study the effects of a collective flow and multiple scattering on two-particle correlation 

o : 

ff^ , measurements in Hanbury-Brown-Twiss intensity interferometry. We find that under a collective 

o 

flow the effective source distribution in a two-particle correlation measurement depends on the 
initial source distribution. In addition, it depends on a collective flow phase function which consists 

Q . 

^ . of terms that tend to cancel each other. As the detected particles traverse from the source point to 



the freeze-out point, they are subject to multiple scattering with medium particles. We examine 



^ ■ the effects of multiple scattering on HBT correlations. By using the Glauber theory of multiple 

scattering at high energies and the optical model at intermediate energies, we find that multiple 
scattering leads to an absorption and an effective density distribution that depends on the initial 
source distribution. 

PACS numbers: 25.75-q 25.75.Gz 25.75.Ld 



I. INTRODUCTION 



Intensity interferometry, as first proposed by Hanbury-Brown-Twiss (HBT) to measure 
the angular diameter of a star using tlie correlation between two photons has been 
applied to optical coherence, subatomic physics, and heavy-ion collisions, utilizing different 
types of particles j^-js^. Much information on the space-time distribution of the emitting 
source is obtained from these measurements. 

Recent experimental measurements of HBT correlations in relativistic heavy-ion collisions 
show only relatively small changes of the extracted longitudinal and transverse radii as a 



function of collision energies 
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The small variation of the HBT radii can 



be explained partly as due to the space-time ordering of the momenta (correlation between 
the position and the momentum) of the produced particles which was used to explain similar 
weak collision-energy dependencies of HBT radii in e~^e~ , pp, pp collisions (see, for example, 
Section 17.4 of Ref. Q). In addition to space-time ordering of momenta, produced particles 
are subject to collective flows. It is important to find out what physical quantities are 
measured by the two-particle intensity interferometry for a chaotic source with a collective 
flow. This will help us understand how collective flows may affect the size parameters 
extracted in two-particle correlation measurements. 

As the detected particles traverse from the source point to the freeze-out point, they 
are subject to final state mean-field interactions and multiple collisions with particles in 
the medium. The effect of final mean-field potential has been studied earlier by Gyulassy 
et al. j?]. They found that for a coherent source the only effect of the mean field potential 
is to redistribute the momentum distribution and the two-particle correlation function is 
not affected. For a chaotic source, the effect of a mean field is given in terms of distorted 
wave functions in the mean field [?[. The effects of the mean field in a chaotic source 
was further studied by Chu et al. |2(| and Shoppa et al. j^]. They found that for low- 
energy detected particles the Coulomb mean-field leads to substantial modification of the 
two-particle correlation, but the effect is small for high-energy detected particles. 

The collisions of the detected particles with medium particles has been considered to be a 
source of chaoticity. It is usually assumed that as a result of the random collisions of medium 
particles, the initial source will evolve into a chaotic source at freeze-out. The source that is 
observed in HBT measurements will be the chaotic freeze-out source, and the HBT radii will 



correspond to those of the freeze-out configuration. For example, in the emission function 
modehng of , an analytical form of the freeze-out phase-space distribution with a 

collective flow is assumed, and the single-particle distribution and two-particle correlations 
are described by a freeze-out distribution. In other examples, hydrodynamics or covariant 
transport theory is followed until freeze-out and the radii of the freeze-out configuration are 
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In such an 



then extracted to compare with experimental HBT radii 
analysis, there is the outstanding puzzle that experimental measurements give -RoutZ-Rside ~ 
0.9 — 1.1 whil e hy drodynamical predictions yield -RoutZ-Rside significantly larger than 1.0 at 



freeze-out l3i 
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The covariant transport theory can explain the magnitude of R gut by 
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assuming a large opacity, but there remains the puzzle that -Rside is under-estimated 

In the context of intensity interferometry, the question whether the freeze-out configura- 
tion in a collective flow is a chaotic source needs to be re-examined carefully as the problem 
of multiple scattering must be treated properly. Because the Hanbury-Brown-Twiss inten- 
sity interferometry is purely a quantum-mechanical phenomenon, the problem of multiple 
scattering must be investigated within a quantum-mechanical framework. We shall focus our 
attention to high- and intermediate-energy detected particles that have sufficient energies 
to propagate from the production point to the detection point. It is necessary to study the 
interference of waves using the probability amplitudes in the multiple scattering process, 
instead of the conventional description of incoherent collisions in terms of probabilities and 
cross sections. The Glauber theory of multiple scattering 4J| has been shown to be a valid 



description for the interaction of a pion with the nuclear medium at a 



300 to 1200 MeV 
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has been found to give a ve rv g ood description of the interaction of a pion with a nucleus 
from 120 MeV to 766 MeV [^. They can be applied here to describe the probability am- 
plitudes for the propagation of energetic detected particles (pions from 0.1 to 1 GeV, say). 
By using the Glauber theory of multiple scattering at high energies and the optical model 
at intermediate energies, we find that the multiple scattering process leads to an absorption 
of particles and an HBT effective density distribution that depends on the initial source 
distribution rather than the freeze-out density distribution. 

Previously, the intensity interference was studied in the framework of probability ampli- 
tudes for the propagation of produced particles We shall used the same framework to 
investigate here the effects of a collective flow and multiple scattering. We examine flrst the 



single-particle distribution in Section II. The probability amplitude for the propagation of 
a produced particle from the production point to the detected point is obtained. The sum 
over the probability amplitudes from all the source points then leads to the single-particle 
distribution which depends on the initial source distribution. In Section III, we study the 
two-particle momentum distribution. For a chaotic source, the quantum effect of symmetriz- 
ing the amplitude for bosons (or anti-symmetrizing the amplitude for fermions) with respect 
to the exchange of the two particles leads to the phenomenon of intensity interferometry. 
We show that under a collective flow the two-particle correlation function depends on the 
initial source distribution and a collective flow phase function. In Section IV, we study the 
effects of multiple scattering of the detected particles with particles in the medium. In Sec- 
tion V, we discuss the relevance of the present results to HBT measurements in high-energy 
heavy-ion collisions. 

II. SINGLE-PARTICLE DISTRIBUTION 

Following arguments similar to those in Ref. j^, we consider the production of a particle 
with four-momentum k = {k,, k^) at a source point x = {x, t) and its subsequent detection 
with four-momentum k = (fc, at the space-time point Xd = {xd,td), as shown in Fig. 
1. For convenience, we shall use the source center-of-mass system as the reference frame to 
measure all momenta and space-time coordinates. 

To describe the single-particle distribution and intensity interference, we need the prob- 
ability amplitude for the particle to be produced at the source point x, to propagate to the 
freeze-out point Xf, and to arrive at the detection point Xd- We can characterize the pro- 
duction probability amplitude for producing a particle of momentum at x by a magnitude 
A{k,x) and a phase (f)o{x). Without loss of generality, the functions A{k,x) and (poix) can 
be taken to be real, and A{k,x) can be taken to be non-negative. The magnitude A{k,x) is 
related to the initial phase-space distribution fi^it{n{x), x) (see Eq. (fT5|) below). The produc- 
tion phase (po{x) describes the degree of coherence or chaoticity of the particle production 
process. We shall be interested in an initially chaotic source which can be represented by 
random and fluctuating production phases (poix) at the source points. 

The complete probability amplitude "^{nx — > kxd) for a particle of momentum k to be 
produced from the source point x, to propagate along the classical trajectory, and to arrive 
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FIG. 1: A particle of momentum k is emitted at a typical source point x of an extended source, is 
subject to a collective flow to reach the freeze-out surface a.t Xf, and is detected with momentum k 
at the detection point Xd- The straight line joining x to Xd is the trajectory of the particle from x to 
Xd- The figure is not drawn to scale. The distance between x and Xd is many orders of magnitude 
greater than the linear dimension of the extended source. 

at Xd with momentum k is 

^!{Kx^kxd)=A{Kx)e'^''^''^{Kx kxd) , (1) 

where iplnx^kxd) is the probability amplitude for the propagation of the produced particle 
from the production point x to the freeze-out point xj and the detection point Xd- From 
the path-integral method, this amplitude is given by [2] 

'iIj{kx — ^ kxd) = ^(classical path, kx kxd) = exp{— z / k{x') ■ dx'}, (2) 

J X 

where the integration is carried out along the classical trajectory, and we have used the 
notation hi{x') to represent the momentum of the particle at the space-time point x' = {x't') 
with t <t' < td- We consider a Lagrangian picture of collective motion in which we follow the 
fluid element at x' and its velocity field I3{x'). For a given detected particle with momentum 
k = {k^, k), the initial momentum k depends on the collective flow velocity f3{xf) at freeze- 
out (see Eqs. Q, Q, and (0) below). The freeze-out coordinate xj in turn depends on the 

initial snnrrp mnrrlinatp t" TVip>rp>fnrp> k.( t"^ is a fnnptinn nf tVif initial sniirpp pnnrrlinatp t" 



The momentum of the produced particle k{x') in different space-time regions along the 
classical path from x to Xd is given in the non-relativistic case by: 

f k(x) + ml3(x'), for t<t' <tf, 

[ k{x) + m/3{xf) — k, for tf <t' < ta, 

and 

r [k{x) + ml3{x')f/2m, for t<t'<tf, 

K^ix') = { (4) 
[ [k{x) + mf3{xf)]'^/2m = k^, for tf<t'< td- 

In the relativistic case, the momentum k{x') under a collective flow with velocity /3(x') is 
given by 



r A(/3(x'))/€(x), for t <t'<t/, 

\ A(/3(a;/))ft;(a;) = k, for tf <t' < ta, 



(5) 



where A(/3(a;')) is a matrix with elements 



Ao(/3) = 7 = 1/Vl - I^P' 

A'o(/3) = A°(/3)=7/3^ (6) 
AU/3) = ^1 + 



and i, /c = 1, 2, 3. 

The momentum of the particle does not change after freezing out at t' — tf. Thus, the 
exponential function in ip{KX kxd) can be separated into the contribution from x to Xf 
and from Xf to Xd, 

%l){nx — > kxd) — exp{— i / n{x') ■ dx' — ik ■ {xd — Xf)}. (7) 

Jx 

In addition the production from the source point x, the particle can also be produced 
from other source points in the extended source. The total amplitude for the detection of a 
particle at Xd is the sum of the probability amplitudes from all source points. After taking 
into account the production probability amplitude Ae**^ at different source points, the total 
probability amplitude for a particle with momentum k to be produced from the extended 
source and to arrive at the detection point Xd is given by 

: {points} ^ ^d) = kxd) 

X 

= x)e^'^°(^) exp{-i P k{x') ■ dx' - ik ■ {xd - Xf)} . (8) 



The single-particle momentum distribution, P{k), which is the probability for a particle 
of momentum k to be produced from the extended source and to arrive at the detection 
point Xd with momentum k, is the absolute square of the total probability amplitude, 

= r <x')-dx' ~tk-{xd-Xf)}\^. (9) 

We expand the righthand side of Eq. Q into terms independent of ^0(3;) and terms con- 
taining (poi^x). We obtain 

X ^'V 

rXf ryf 

X exp{-z / K{x')-dx' + i K{y')dy'} . (10) 

J X Jy 

For an initially chaotic source, the phases at different source points can be described by a 
random and fluctuating phase (t)o{x) with a period much shorter than the mean-life of the 
source. For such a source, the second term on the righthand side of the above equation gives 
a zero contribution because the large number of terms with slowly varying magnitudes, but 



rapidly fluctuating random phases, cance 



coordinate (see Supplement 17.1 of Ref. 



out one another in the sum over the source time 
^). The properties and the validity of using the 
time average for a chaotic source have been discussed in detail in Chapter X of the text of 

n 

Born and Wolf [531. Therefore, Eq. ^ becomes 

P{k) = YA\K{x),x). (11) 

X 

The summation over the source points necessitates the specification of the initial density 
p{x) of the source points per unit space-time volume at the point x. With this specification, 
the summation should be transcribed as an integral over x, 



Y- ^ d'xpix)... (12) 

X '' 

Therefore, we can rewrite Eq. pi|) as 

P{k) = [ d'^xp{x)A'^{K{x),x) , (13) 



which is independent of Xd- We can compare this equation with the properties of the phase 

snanp rlistriKiitinn fimr'tinn Onp> nan dpsnrWie' flip nKsp>rvp>rl nnrtir-lfis as rrimincr frnm tVip> 



initial phase space distribution fimt{i^{x) , x) defined by 



P{k) = J d'^x finitiK{x),x) 



(14) 



A comparison of the above equation with Eq. (jl3|) shows that the observed particle with 
momentum k comes from the initial phase space distribution 



/init('^(a;),a;) = p{x)A'^{k{x),x) 



(15) 



One can alternatively describe the observed particles as coming from the freeze-out phase 
space distribution ff{k,Xf) defined by 

(16) 



P(k) = j d Xf ffik,Xf) = J d X /init(/«(a;), x). 

The initial distribution and the freeze-out distribution are equivalent descriptions and they 
are related to each other by 



d{xf,xV^ 



fmit{l^{x),x) 



d{x, x°) 



d{xf,xj) 



p{x) (k{x) , x) 



(17) 



where 



d{x,x^)/d{xf,x^j) is the Jacobian determinant arising from the mapping of the 



initial source point {a;,x°} to the freeze-out point {xf,x'l} due to the collective fiow. 



III. TWO-PARTICLE DISTRIBUTION 

We consider the case in which a particle of momentum k,i{xi) starts from xi, propagates 
to the freeze-out point X/i, and arrives at the detection point Xdi with momentum ki, and 
another identical particle of momentum K2{x2) starts from X2, propagates to the freeze-out 
point Xf2, and arrives at Xd2 with momentum k2, as indicated by the solid lines in Fig. 2. 
The probability amplitude for the production of the particle with momentum Kj{xi) at Xi 
is given by A{Kj{xi), Xi)e'^^°^^^\ Therefore, the probability amplitude for the two particles 
to be produced at the source points, to propagate from the source points to the freeze-out 
points, and to arrive at the detection points is 

A(ft:i(xi),Xi)e'<^«("^)A(/€2(x2),X2)e''^°("2) exp{-^ /^/^ /€i(x') ■ rfx' - ^fci ■ (x^.-Xfi)} 

X exp{-i K2{x') ■ dx' - ik2 ■ {xd2 - x/2)} • (18) 



source 




FIG. 2: A particle of momentum k\ is detected at x^i and another identical particle with momen- 
tum k2 is detected at the space-time point Xd2- They are emitted from the source point xi and 
X2 of an extended source. The solid lines joining xi — > x^i — ^ x^i for ki and X2 ^ Xf2 ^ x^s for 
A;2, and the dashed lines joining xi x'j^ — > X(i2 for /c2 and X2 x'^^ ~^ ^di for ki are possible 
trajectories. 

However, this is not the only probability amplitude contribution for two identical particles 
produced from Xi and X2 to arrive at and 2:^2 • The particle of momentum ki detected 
at Xdi can also be produced at X2 with momentum Ki{x2) , propagate from X2 to the freeze- 
out point 2:^2(^2), and arrives at xai with momentum ki, while the other identical particle 
of momentum ^2(2^1) starts at Xi, propagates from xi to the freeze-out point x'fi{xi), and 
arrives at Xd2 with momentum k2, as indicated by the dashed lines in Fig. 2. As the distances 
between the freeze-out points and the source points are much smaller than the distances 
between the freeze-out points, and the detection points and two-particle correlations occur 
for ki close to k2, it is reasonable to make the approximations x'^^ = ^fi and Xf2 = x'^2 (^^^ 
Fig. 2). The probability amplitude for this occurrence along the trajectories x\ x'fi — > Xd2 
for k2 and X2 — > x'f2 — x^i for ki is 

yl(fi:i(a;2),a;2)e*'^°("2)A(fi:2(a;i),a;i)e'^°(^i) exp{-?/;f k^{x') ■ dx' - iki ■ {xai - Xf2)} 

X exp{-i J^/' K2{x') ■ dx' - ik2 ■ {xd2 - x/i)} . (19) 

Rppaiisp of tVip indistincriiisViaKilitv of tVip nartiplps and tVip Rnsp-P,instpin statistips of idpn- 



tical bosons (or the Fermi-Dirac statistics for fermions), the probabihty amphtude must 
be symmetrical (or antisymmetrical) with respect to the interchange of the labels of the 
particles which distinguish them. In this case, the only labels which distinguish the two 
identical particles are the source point coordinates xi and X2, because one particle of mo- 
mentum ki has been determined to have been detected at Xdi and the other identical particle 
of momentum k2 at Xd2- The probability amplitude must be symmetrical for bosons and 
anti-symmetrical for fermions, with respect to the interchange of the labels Xi and X2- Ac- 
cordingly, the probability amplitude which satisfies this symmetry is the sum (or difference) 
of Eqs. (UHl) and (HH) divided by ^2: 

^( A(/€i(xi),Xi)e^'^°(^iM(K2(x2),X2)e*^o(^2)exp{-i r\i{x') ■ dx' - ik^ ■ {xdi-Xfi)} 

V2 I J XT, 

rxf2 

X exp{-i / K,2{x') ■ dx' - ik2 ■ {xd2 - 3^/2)} 
+M(/ti(x2),X2)e^'^"(^2M(K2(xi),a;i)e*'^o(^'i) exp{-i ki{x') ■ dx - iki ■ {xdi - x/2)} 

J X2 

rXfi 

X exp{— 2 / K,2{x') ■ dx' — ik2 ■ {Xd2 — Xfi)}} 

= e*'^°("^)e^'^«("^)$(/€i«:2, X1X2 ^ XdiXd2) , (20) 

where ^ is 1 for bosons, —1 for fermions, and ^{k,ik,2,xiX2 — > XdiXd2) is the part of the 
probability amplitude in the above equation which does not depend on 0o- It is defined by 

$(/«iK2 : X1X2 XdlXd2) 

= —i=<.A{Ki{xi),Xi)A{K2{x2)x2)eyip{—i / k,i{x') ■ dx' — iki ■ {xdi — Xfi)} 
v2 I Jxi 

rxf2 

X exp{-'i / K2{x') ■ dx' - ik2 ■ {xd2 - 3^/2)} 

J X2 

rxf2 

+ 6A{ki{x2), X2)A{k2{xi), Xi) exp{-z / ki{x') ■ dx' - iki ■ {xdi - Xf2)} 

Jx2 
rXfi 

X exp{-i / K2{x') ■ dx' - ik2 ■ {xd2 - Xfi)} . (21) 

J Xl 

Besides originating from the source points Xi and X2, the two particles can also be produced 
at other source points in the extended source. The total amplitude is the sum of amplitudes 
from all combinations of two source points. Therefore, the total probability amplitude for 
two identical particles to be produced from two source points in the extended source and to 
arrive at their respective detection points Xdi and Xd2 with momenta ki and /c2 is 

^(fcifc2:{pointr}^a;,iXd2)= E e^'^^^^'^^e^'^''^^^) [ ^KiK2:xiX2-^XdiXd2) 

{xi,X2} 



where the sum is carried over distinct combinations of {xi,X2}- The two-particle momentum 
distribution P{ki, ^2) is defined as the probabihty distribution for two particles of momenta 
ki and k2 to be produced from the extended source and to arrive at their respective detection 
points Xdi and Xd2- From Eq. (j^ . it is given by 

P{k,,k2) = ^|^K«:2:{Sirn (23) 

We now apply the results of Eqs. ()21|) - (j^n|) to examine the momentum correlations of a 
chaotic source. A chaotic source is described by a phase function (f)o{x) that is a random 
and fiuctuating function of the source point coordinate x. For a chaotic source, we again 
make use of the random and fiuctuating nature of the phases by substituting Eq. 1)2111 into 
Eq. ()22j) and expand the righthand side of Eq. ()23p . We separate out terms which are 
independent of (po and terms which contain 0o- We obtain 

P{ki,k2) = ^ V($*(/tiK2:yil/2^a;diXd2)|w=-i$(KiK2:xiX2->XrfiXrf2) 

+ {K1K2 : yiy2-^XdlXd2) I !'2=-i ^{kiK2 : XiX2^XdlXd2jf 

yi=x2 ') 



+ 2 ^ 

{xix2}7^{yiy2} 



X <I>*(KiK2 : yiy2-^XdlXd2)^{'ilti2 ■ XiX2^XdlXd2)^ ■ (24) 



The two terms in the first summation on the righthand side are equal because of exchange 
symmetry. For the chaotic source, the last term in Eq. gives a zero sum because the 
contributions of a large number of terms with similar magnitudes but random and fiuctuating 
phases cancel out. Therefore, Eq. becomes 

P{ki, k2) = |$(kiK2 : X1X2 XdiXd2)\'^ , (25) 

Xl,X2 

Converting the summations in Eq. into integrals with the transcription p2j) . we can 
rewrite the total probability as a double integral over the source point coordinates Xi and 

X2, 

P(k„k,) = p(x,)pfe)H.(.,.. : x.x. ^ (26) 



In the above equation, |$p can be obtained by using Eq. (j2H). We get 

|$|^ = A\ki{xi),Xi)A'^{k2{x2)x2) 

+6A{ki{xi),Xi)A{k2{x2)x2)A{ki{x2),X2)A{k2{xi),Xi) , (27) 



where 



i-Xfi j-xf2 rxfi j-xf2 

A = — j Ki{x') ■ dx' — / K2{x') ■ dx' + / K2{x') ■ dx' + / Ki{x') ■ dx' 

J Xi J X2 Jxi J X2 

^\ki{x') - K2{x')] ■ dx' - n\K2{x') - Ki{x')] ■ dx' , (28) 

Xi J X2 



and 



B = ki- (X/i - Xdl) + k2 ■ {Xf2 - Xd2) - ki ■ (x/2 - Xdl) - k2 ■ {Xfi - Xd2) 

= {ki- k2) ■ Xfi + {k2- ki) ■ Xf2. (29) 



Therefore, we have 



P{ki,k2) = J d^Xip{xi)A'^{Ki{xi),Xi) J d'^X2p{x2)A'^{K2{x2),X2) 

+e I rf^xip(xi)A(Ki(a;i),a;i)A(/«2(xi),xi)e'('=i-'=2)-"i exp{-t £'\[ki{x') - K2ix')] - [ki - A^a]} ■ dx'} 

X / d'x2p{x2)A{K2{x2),X2)A{K,{x2),X2)e-'^'''-'''^-'''exp{-l r\[K,{x') - K2{x')] - [k, - k2]} ■ dx'} 
J JX2 

(30) 

Using Eq. (jTHl), we can rewrite the above equation as 

P{ki,k2)=P{ki)P{k2) 



+e 



^4^g*{fci-fe2) xg^p|_^ P{[ki(x') - K2{x')] - [ki - k2]} ■ dx'}p{x)A{Ki{x),x)A{K2{x),x) 



(31) 



which is independent of Xd- It is convenient to rewrite the two-particle distribution function 

as 

P{ki, k2) = P{ki)P{k2) (l + ^ / e'('=l-^2)•x+^0.{.,A,fc2)^^^(^. ^^^2) ^ (32) 
where pes is the effective density 

p{x)A{ki{x),x)A{k2{x),x) 
Pes{x; kik2) = , , (33) 

JP(k,)P(ko) 



and (pc{x,kik2) is the collective flow phase function 

kik2) = - P{[Ki(a;0 - K2{x')] - [ki -k2]}- dx' . (34) 

The HBT two-particle correlation function R{ki, /C2) = P{ki, k2) / P{ki)P{k2) — 1 becomes 

R{ki, k2)=e I d^a;e*('=i-'=2>^+*^=(^''=^^2Veff(a;; /ti^a) \ (35) 

From Eq. 0, we have 

[K^{x') - K2{x')] - [k^ - k2] = [A(/3(x')) - AH3ixf))]A-\f3ixfmk^ - k2). (36) 

The collective flow phase (pdx-, ^^1^2) becomes 

U^Mk2) = - r dx' ■ {[Ai(3ix')) - A{(3{xf))]A-\(3ixj)){k, - k2)} . (37) 

Thus, under a collective flow the two-particle correlation function depends on the initial 
source distribution and a collective flow phase function (p^x, kik2), as given by Eqs. (jH^ - 
When there is no collective flow, 4>c{x, /C1/C2) is zero, and P{ki, k2) and -R(fci, fco) reduce 
to the usual two-particle correlation functions of Eqs. (17.29) and (17.38) in Ref. 2|. 

We note that the collective flow phase function 0c (a;, kik2) consists of the differences 
A{(3{x')) — A{f3{xf)) and ki — k2 which can lead to substantial cancellation. In fact, if l3{x') 
reaches the asymptotic velocity f3{xf), then there will be no contribution from this path 
element dx' at x' to the phase function (pdx, kik2). 

The velocity proflle f3{x') of a source can be obtained by solving the relativistic hydrody- 
namical equations. We can get some ideas on the expansion scenario for a highly compressed 
relativistic ideal fluid initial at rest from the results of Rischke and Gyulassy 3^ . We can 
consider a spherical flreball in which the collective flow is along the radial direction. (The 

n 

transverse expansion of a Bjorken cylinder can be discussed in a similar way [23|-) When 
this highly compressed fluid is allowed to evolve hydrodynamically, fluid elements at the 
surface reach their asymptotic velocity (close to the speed of light) very rapidly — nearly 
at the beginning of the expansion. Fluid elements at the central region remain essentially 
at rest until they are carried outward by the strong rebounding wave, and they reach the 
asymptotic velocity also rapidly. 

From Fig. 2 of the velocity of the fluid element initial at (r, 9, (p) is approximately 



where tc = t + 2{R — r), r < R, and R is the radius of the initial source. The trajectory for 
the fluid element is r'{t') = r + i3{x')t', 0'{t') = 6, and (j)'{t') = cf). For such a collective flow, 
we can describe this as an initial source up to t = tc, and the source undergoes collective 
flow at t > tc with a flow velocity reaching the asymptotic velocity approximately in a 
step-wise manner. Fluid elements reaching the asymptotic velocity do not contribute to the 
collective phase (^d^, kik2), as one notes previously in Eq. (jHTj) . The collective flow phase 
function (pd^, kik2) is approximately zero for this case. For other initial conditions and 
compression, the radial expansion and rapid rise of the collective velocity to the asymptotic 
velocity should be similar, as the collective motion arises from the expansion of a compressed 
matter into an empty vacuum. The magnitude of the asymptotic velocity may depend on 
the degree of compression (and therefore on the heavy-ion collision energy). The collective 
flow phase function 0c will likely remain small because of the substantial cancellation in 
Eq. (j37|) mentioned above. We should expect that the effective density measured in HBT 
measurements should depend essentially on the initial source distribution. Much more work 
will be required to examine further the behavior of (pd^, ^1^2) and R{kik2) for general flows 
and initial conditions. 



IV. EFFECTS OF MULTIPLE SCATTERING AND OPTICAL POTENTIAL 



The observation of the intensity interferometry depends on the degree of coherence or 
chaoticity of the source. As the initial source particles traverse from the source point Xi 
to the freeze-out point xji under a collective flow, they are subject to multiple scattering 
with medium particles. One may think that as a result of these scatterings, the source 
becomes chaotic at freeze-out and the distribution observed in an HBT measurement should 
correspond to the freeze-out distribution 

The problem of multiple scattering must however be treated properly in the context 
of intensity interferometry. As the Hanbury-Brown-Twiss effect is purely a quantum- 
mechanical effect, the problem of multiple scattering must be investigated within the 
quantum-mechanical framework in terms of probability amplitudes instead of incoherent 
collisions involving probabilities and cross sections. 

Within such a theoretical framework, we can study our problem in the Glauber model 



4J, |45[ or the optical model 
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depending on the energy of the detected 



particle relative to the medium. The Glauber theory has been applied successfully to study 
the interaction of a pion with the nuclear medium at a pion energy from 300 to 1200 MeV 



45 



an therefore be used to study the interaction of an energetic pion with 
the hadron medium produced in high-energy heavy-ion collisions. At these energies, the 
trajectory can be approximately described as a straight line, the probability amplitude for 
the particle to traverse from x to x' after suffering multiple scatterings with particles in the 
medium is 

N 

exp{z</)^(a; x')} = exp{z ^ x{x'^ - Sji_)e{x\ - Sj\{)e{sj\\ - x{)], (39) 

where the function ~ ■Sjx) is the Glauber phase shift function for the interaction of the 
produced particle with a medium particle at Sj and the sum is carried over all medium 
particles having coordinates Xj = {tj, Sj\\, Sjj_),j = 2, 3, [44j. The phase shift function 
can be obtained from two-body scattering data and is an analytical function of the 
transverse coordinate b. The subscript _L in Eq. denotes the component transverse to 
the particle trajectory and the subscript || denotes the component along the trajectory. 

The wave function Eq. (j39p from Glauber multiple scattering theory contains a wealth 
of relevant information. It depends on the coordinates of all the particles with which the 
incident particle has interacted. It treats correctly the case of no scattering and multiple 
scattering, even up to the extreme case of — 1 scatterings in succession. It makes no 
difference whether the medium particles are dense in close proximity or dilute in far separa- 
tion. Information on the density of medium particles can be provided when one integrates 
out the distribution of the medium particles. 

One can show that starting with the multiple scattering wave function Eq. ()39|) . one can 
construct the Wigner function in the transverse and longitudinal degrees of freedom. The 
multiple scattering wave function Eq. (j39|l leads naturally to to the diffractive transport 
theory, and the equation of motion for the Wigner function is just the Boltzmann equation 
[5^ . The solution of the Wigner function can be decomposed into multiple-scattering com- 
ponents each of which corresponds to the scattering of the incident particle with a definite 
number of medium particles. The transverse width of the n-scattering component increases 
with ^/n and the maximum of the longitudinal distribution of the n-scattering component 
is located longitudinally at n times the mean-free path. 

As f^laiiKfr mnltir>l(=> spatt.princr tVipnrv crivps tVi(=> plassiral HifFraptivp transnnrt tVipnrv nn 



the one hand, and retains the wave nature of the propagation on the other hand, it is 
therefore appropriate to use Glauber multiple scattering wave function to investigate the 
effects of multiple scattering in intensity interferomctry. 

Using this description for the multiple scattering of the particle, the probability amplitude 
for the produced particle to travel from the source point x to xa becomes 

'^{Kx^kxd)=A{Kx)e'-'^'''^^ip{KX kxd) (40) 

where (f){x) is 

= (j)o{x) + (j)s{x Xf), (41) 

00 (2^) is the initial random and fluctuating phase due to the chaoticity of the production 
source, and 0s(x Xf) is 

N 

<t>s{x Xf) ^Y.x{xf± - Sj^_)e{xf\\ - Sj{)e{Sj\^ - x{). (42) 

As Xf is a function of x, the phase function (j)s{x Xf) is a function of x and can be 
abbreviated as 0s (x). It contains a real and imaginary part. The imaginary part is given by 

N 

Im (f)s{x) = ^Xm x{xf± - Sj±)9{xfii - Sj||)^(Sj|| - ajy), (43) 

and it represents the degree of absorption as the detected particle passes through the 
medium. 

With this modification, the derivation of the single-particle and two-particle distributions 
can be carried out as in the last sections. As the multiple scattering bring in a complex 
phase (psix) that depends on x, the quantum description of the scattering amplitude leads 
to the single-particle distribution 

P{k) = j dX2dX3dX4...dXN J2 e'^'*''^''^~''''^^''^^A'^{K{x),x)p^ed{x2: X3, X4, Xn) , (44) 

where Pjned{x2, X3, X4, Xn) is the distribution of the medium particles. Note that the real 
part of (f)s{x) cancel out in the factor 0s (x) — 0s (^) above equation. Consequently, 

when the multiple scattering between the detected particle and the medium is taken into 
account, the single-particle distribution for a collective flow is 

P(k) = I d^x ~^'^''^p{x)A^U{x),x) , (45) 



where 



,X3,X4, ...,Xn) 



In this case, we have 



P{k) = / d^a;/init(K(x),x), 



where 



/init(/t(a;),x) = e' 



-2 Im 4>s{x) ( \ a2 



p{x)A {k,{x), x) 



(46) 



(47) 



(48) 



The two-particle distribution for a system with a collective flow can be obtained as in the 
last section, and we get 



P{h,k2) = P{ki)P{k2) 



+ e 



^4^gi(fci-fc2)-x+i<^c(a;,A:ifc2)g-2 Xm U^) p(^x) A{Ki{x) , x) A{k2{x) , x) 



or 



p{k^,k2) = P{k^)P{k2)[i + e 



^i{ki-k2)-x+i<j)c{x,kik2)^-2 Xm <f>s{x) 



pesix; kik2) 



(49) 



(50) 



where 0c(a^, ^1^2) is given by Eq. (jH^ and Pefr(a;; ^1^2) is given by Eq. Thus, multiple 

scattering leads to an absorption of produced particles as they traverse from the source point 
to the freeze-out point. There is an absorption factor in both the single-particle distribution 
and the two-particle correlation function. The effective density distribution pefi revealed by 
HBT two-particle correlation measurements depends on the initial source distribution and 
not necessarily on the freeze-out density distribution. 

At lower energies, the interaction of the pion with the medium can be described by an 

42, H H H H Q- 



optical model 



For example, in the interaction of a pion with a 



nucleus, a local phenomenological optical potential has been successfully applied to explain 
the TT- nucleus elastic scattering data from 120 MeV to 766 MeV |4^ |42l |^ . The effects of 
intermediate resonances such as A(1232) shows up as giving rise to a peak in the imaginary 
part of the optical potential near the A(1323) resonance energy. The optical potential is 
normally non-local but can be described in terms of a local-equivalent potential by redefining 
some kinematic quantities Q]. A "model-exact' microscopic description of the optical po- 



in terms of the interaction between the pion and medium particles, including the effects of 



HQ- 



the A(1323), Ai3(1520), and Fi5(1680) resonances 

Based on the successes of the optical model in describing the interaction between a pion 
and nuclear matter, it is reasonable to use similar descriptions to study the dynamics of an 
intermediate-energy pion in hadronic matter produced in high-energy heavy-ion collisions. 
An optical model is specially appropriate as the Hanbury-Brown-Twiss intensity interfer- 
ometry arises from an optical interference of the wave amplitudes. 

One can accordingly introduce a complex optical potential V{k, x) with a negative imagi- 
nary part to describe the interaction between the detected particle and the medium particles 
as it travels from the source point to the freeze-out point. The optical potential V depends 
on the momentum of the particle k. The momentum dependence is particularly pronounced 
in the neighborhood of a resonance. In principle, the optical potential can be determined 
from two-body data and the distribution of particles in the medium, as carried out for exam- 
pie . Ref. U.*. tKe .t.acUon of a. op*a. po.entia. the a.p...de 

for a produced particle to propagate from x to x" is given by [44 1 

rx" 1 

exp{i(f)s{x ^ x")} = exp{—i / -V{k,x')dx'u}, (51) 

Jx V 

where v is the magnitude of the velocity of the detected particle relative to the hadron 
medium along the direction of propagation. Therefore, when the multiple scattering between 
X and X/ is taken into account, the probability amplitude for the particle to propagate from 
X to Xd becomes 

^{Kx^kxd)=A{Kx)e'^'-''^-^{Kx kxd), (52) 

where 

</,(x) = (f)o{x) - r" -V{k, x')dx\. (53) 

Jx V 

The phase shift due to the optical potential, (f)s{k,x), is a complex quantity containing a 
real and an imaginary part, 

(l),{k,x) = - J^^ ^V{k,x')dx\^. (54) 
As Xf is a function of x, the phase function 0s(x — > Xf) is a function of k and x and can be 

aKKrp>viatp>rl Ptfi rhJh. t\ 



With this modification, the derivation of the single-particle and two-particle distributions 
can be carried out as in the case with the Glauber theory. The two-particle correlation 
function R{ki,k2) becomes 

2 



(55) 



where 

pxd / I 1 \ 

Tie Usiki, x) - 0*(A;2, x)) = - —Tie V(ki, x')dx\, 7^e V(k2, x')dx'u , (56) 

Jx \Vi " V2 V 

Im {(j)s{ki^x) — (j)*{k2^x)) = — I — Xm V{ki,x')dx'u H Tm V{k2,x')dx'u ) 

= 2Im (j)s{x) (57) 

and Vi = \{ki)\\\/k^. For low-energy detected particles, the velocities of the two detected 
particles can differ substantially and the real part of the phase function does not cancel 
exactly. There will be a substantial correction to the two-particle correlation function at 
low energies 



pions particles with kinetic energies close to or greater 
than its rest mass, in which we shall focus our attention, fj ~ 1. The real parts of the 
phase function (f)s{k,x) cancel and there will be negligible modification of the two-particle 
correlation function. We again obtain the results of Eqs. (jUj), (jlH)), and (^Hj) for the single- 
particle distribution, and the result of Eq. for the two-particle distribution. 

In applying the Glauber theory at high energies or the optical model at intermediate 
energies, the precise energy at which one needs to switch from the Glauber theory to the 
optical model needs not concern us here at present as both descriptions lead to a complex 
phase shift function (psix). The imaginary part of the phase shift function gives rise to an 
absorption. The real part of the phase shift leads to negligible modification of the two- 
particle correlation function. The conclusion concerning the real phase of coordinate is a 
rather general result. Any final-state interaction, that leads to an additional real phase 
function of coordinate and a weak dependence on momentum, does not modify the two- 
particle correlation. As a consequence, the effective density distribution depends on the 
initial source distribution and not necessarily on the freeze-out distribution. The intensity 
interferometry arises from the difference in the phases when the two particles takes on two 
different sets of histories. For either set of histories depicted in Fig. 2, even though the 

nmr>li+iirlfi fnr thp r>rnr>n crnfirm frnm +n Tj'i rlp>r>pinrls nn +Hp> r>nsifirms nf thp mprliiiTn 



scatterer which can be random, the real parts of the phases in propagating from the source 
point to the freeze-out points are the same. They cancel out when we evaluate the HBT 
two-particle correlation function R{ki, k2). 

Previously, Gyulassy et al. found that the two-particle correlation function for a coherent 
source is unaffected by a final-state mean-field potential. For a chaotic source, which is the 
object of our interest here, they expressed the two-particle correlation function in terms 
of distorted waves in the mean-field [Eq. (5.32) of j7|]. The present results for the optical 
potential corresponds to an explicit evaluation of this two-particle correlation function using 
the eikonal approximation for high-energy particles. In the process of this evaluation, we 
find that the phase distortion due to a real mean-field potential cancels out, as we remarked 
earlier, and lead to negligible effect on the two-particle correlation at high particle ener gies . 



The distortions can however be substantial for low-energy particles as discussed in 
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V. CONCLUSIONS AND DISCUSSIONS 



We start from the probability amplitude written in terms of the path integral over the 
classical trajectory for a single particle and a pair of identical particles. The sum over the 
probability amplitudes for a chaotic source leads to an intensity interference for the detection 
of two identical particles. 

As the detected particles traverse from the source point to the freeze-out point, they 
are subject to scattering with particles in the medium. We have examined the effects of 
multiple scattering. Because the Hanbury-Brown-Twiss intensity interferometry is purely a 
quantum-mechanical phenomenon, we investigate the problem of multiple scattering within 
the quantum-mechanical framework of the Glauber theory and the optical model. We find 
that multiple scattering leads to an effective density distribution that depends on the initial 
source distribution rather than the freeze-out density distribution. This conclusion follows 
from the wave nature of the detected particle, which is an important ingredient for the occur- 
rence of intensity interferometry. The Glauber multiple scattering amplitude and the optical 
model scattering also represent coherent scattering processes j3] as the corresponding phase 
function is an analytical function of the coordinate. There is the expected absorption arising 
from the imaginary part of the phase shift function in the multiple scattering process. 

We have studied the multiple scattering process at high energies using the Glauber model 



and at intermediate energies using the optical model in the eikonal approximation. At low 
energies when these simple approximations may not be valid, it will be of great interest to 
study a quantum mechanical description of the multiple scattering process using either a 
quantum mechanical many-particle theory or a pion mean-field optical potential without the 
eikonal approximation. 

In the presence of a collective flow, we find that two-particle correlation measurements 
lead to an effective source distribution that depends on the initial source distribution. The 
relevant momentum is the initial momentum shifted from the detected momentum downward 
by the collective flow, as given by Eq. (jHHj) obtained here. In addition, the collective flow 
leads to a phase function (pd^, kik2) in Eq. ()34|) which will modify the effective two-particle 
correlation. It consists of terms which tend to cancel each other. We examine sample 

n 

results of hydrodynamical calculations of Rischke and Gyulassy [33|| for a highly compressed 
relativistic fluid initially at rest. We find that the flow velocity of a fluid element reaches 
the asymptotic velocity rapidly at the surface, or in the interior when the rebounding wave 
carries the fluid element outward. There can be a substantial cancellation in the collective 
flow phase function (pdx, kik2)- As a consequence, the effective density measured in HBT 
measurements is expected to depend essentially on the initial source distribution for such an 
expansion. Much more work will be required to examine further the behavior of (pdx, kik2) 
and R{ki, /C2) for general flows and initial conditions. 

While further studies are continuing, it is interesting to explore the possibility that HBT 
measurements are indeed related mainly with the distribution of the initial source distribu- 
tion. In that case, we expect that the initial source transverse dimension should be given 
approximately by the spatial dimension of the colliding nuclei which should be nearly inde- 
pendent of the collision energy. We also expect that in the initial source prior to the collective 
flow, the transverse size in the "out" direction should be approximately the same as that 
in the "side" direction. Hence, the HBT transverse radii should be approximately the same 
as the colliding nuclear radii, have only a weak collision-energy dependence, and Ront/Rside 
should be approximately close to 1. These expectations are consistent with the gros s features 



of HBT transverse radii in high-energy heavy-ion collisions 
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